import numpy as np
project_name=__name__.split('.')[-1]
# 时间数组
t_array_target = np.linspace(0.0, 20.0, 30)  # 从0到20秒，生成30个点


# 正弦函数参数
p_min = 10e6     # 压力最小值
p_max = 20e6     # 压力最大值
phase_start = 0          # 起始相位（弧度）
phase_end = 2 * np.pi    # 结束相位（360度）

# 将时间数组映射到 [phase_start, phase_end] 的相位区间
phase_array = np.linspace(phase_start, phase_end, len(t_array_target))

# 使用正弦函数生成压力曲线：从 0° 到 360°，值域从 0 到 1
# 这里使用 sin(x) 在 [0, 2π] 区间的变化趋势
p_array_target = p_min + (p_max - p_min) * (np.sin(phase_array) + 1) / 2

t_array_target=np.append(t_array_target,21.0)
p_array_target=np.append(p_array_target,0.0)

cstar = 1580
rho = 1780
n = 0.35
At = 0.01